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Abstract 

In a system of atoms with large positive scattering length, weakly-bound diatomic molecules 
(dimers) are generated dynamically by the strong interactions between the atoms. If the atoms 
are modeled by a quantum field theory with an atom field only, condensates of dimers cannot be 
described by the mean-field approximation because there is no field associated with the dimers. 
We develop a method for describing dimer condensates in such a model based on the one-particle- 
irreducible (1PI) effective action. We construct an equivalent 1PI effective action that depends not 
only on the classical atom field but also on a classical dimer field. The method is illustrated by 
applying it to the many-body behavior of bosonic atoms with large scattering length at zero tem- 
perature using an approximation in which the 2-atom amplitude is treated exactly but irreducible 
iV-atom amplitudes for N > 3 are neglected. The two 1PI effective actions give identical results 
for the atom superfluid phase, but the one with a classical dimer field is much more convenient for 
describing the dimer superfluid phase. The results are also compared with previous work on the 
Bose gas near a Feshbach resonance. 

PACS numbers: 03.75.Nt, 34.50.-s, 21.45. +v 
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I. INTRODUCTION 



The development of the technology for cooling atoms to ultralow temperatures led to 
the discovery of Bose-Einstein condensation in dilute atomic gases P, 0, 3 . It opened up 
the possibility for experimental study of Bose gases. By controlling the scattering length 
of the atoms using Feshbach resonances, it is possible to study strongly-interacting Bose 
gases 0,0,0]. Among the interesting phenomena that have been observed is atom- molecule 
coherence 0, 0, 0] , which is characterized by oscillations whose frequency coincides with the 
binding energy of a diatomic molecule. 

There have been several studies of a Bose gas of atoms near a Feshbach resonance. Radz- 
ihovsky, Park, and Weichman [10;] and Romans, Duine, Sachdev, and Stoof [ll[ predicted 
that in addition to the normal phase and an atomic superfluid phase, there is a molecular su- 
perffuid phase. In the atomic superfluid phase, there is a Bose-Einstein condensate of atoms 
and there may also be a Bose-Einstein condensate of diatomic molecules (dimer condensate). 
In the molecular superfluid phase, there is a dimer condensate but no atom condensate. Lee 
and Lee used the renormalization group to study the nature of the phase transition between 
the atomic superfluid and the molecular superfluid Basu and Mueller studied various 
instabilities of the superfluid phases 13]. The studies of Refs. |To[ and 11] indicate that 



at zero temperature, there is a quantum phase transition between the atomic superfluid 
and the molecular superfluid. The relevant interaction parameter is the detuning energy 
of the molecule or, equivalently, the scattering length a of the atoms. The quantum phase 
transition occurs at a positive critical value a c of the scattering length that depends on the 
number density. 

A Feshbach resonance allows the scattering length a to be tuned arbitrarily large. When 
the scattering length a is large compared to the range (and all other length scales set by 
the interactions), identical bosons have universal few-body properties that are determined 
by a and governed by a discrete scaling symmetry [3]. Efimov showed in 1970 [15] that in 
the unitary limit a = ±oo, there are infinitely many arbitrarily shallow 3-body bound states 
(called Efimov trimers). Their binding energies have a geometric spectrum: 

4 n) = 22.7- 2{n - n ^h 2 K 2 Jm , (1) 

where m is the mass of the bosons and is the binding wavenumber of the Efimov trimer 
labelled by n*. We will refer to /t* as the Efimov parameter. The spectrum in Eq. ([T]) is 
consistent with a discrete scaling symmetry with discrete scaling factor e 7 ^ ~ 22.7, where 
Sq ~ 1.00624 is a transcendental number. Efimov showed that the discrete scaling symmetry 



is also relevant when a is large but not infinite [161 . |17| ] . For example, few-body observables 
involving only particles with zero momentum must have the form a p f(so\og(aK*)), where 
p is the power required by dimensional analysis and f(8) is a dimensionless function with 
period 2ir. We will refer to such log-periodic behavior as Efimov physics. 

In the previous analyses of the strongly-interacting Bose gas in Refs. 0, H E O 



Efimov physics has been completely ignored. An obvious qualitative consequence of Efimov 
physics is the possible existence of a trimer superfluid phase in which there is a Bose- 
Einstein condensate of Efimov trimers (trimer condensate), but no atom condensate or dimer 
condensate. The methods used in Refs. [l(J 11, [H, EH would be incapable of describing 



a trimer superfluid phase. In these analyses, the model used as a starting point was a 
quantum field theory with quantum fields i[>a{t, t) and i^m{^, t) that annihilate an atom and 



a molecule, respectively. Most of the quantitative analyses in Refs. [10|, [TJJ were carried 
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out using the mean-field approximation. The order parameters associated with the atomic 
superfluid and molecular superfluid phases were the mean fields (V>a) and (V>m), respectively. 
In the atomic superfluid phase, (V>a) 7^ and (i^m) 7^ 0, while in the molecular superfluid 
phase, (iI>m) 7^ and (tpA) = 0. If the scattering length in this model is tuned to be 
sufficiently large, Efimov trimers will be generated dynamically by the strong interactions 
between the atoms. However the mean-field approximation cannot reveal the existence of 
a trimer superfluid phase, because there is no quantum field whose expectation value can 
indicate the presence of a condensate of Efimov trimers. 

There is an analogous problem with the dimer superfluid phase in a quantum field theory 
with only an atom field ^(r, t). If the scattering length a is positive and sufficiently large, 
the dimer is generated dynamically by the strong interactions between the atoms. The mean 
field {ipA} provides an order parameter for the atom superfluid phase. However the mean- 
field approximation cannot reveal the existence of a dimer superfluid phase, because there 
is no quantum field whose expectation value can indicate the presence of a condensate of 
dimers. 

In this paper, we show how a dimer superfluid phase can be studied in a quantum field 
theory in which the only quantum field is the atom field ^(r, t). Our basic tool is the 
one-particle-irreducible (1PI) effective action T[A], which is a functional of a classical atom 
field A(r,t). This effective action encodes complete information about the quantum system. 
The variational equations for the 1PI effective action provide a straightforward description 
of states containing an atom condensate or a mixture of an atom condensate and a dimer 
condensate. They also contain complete information about states containing a dimer con- 
densate but no atom condensate, but this information is encoded in a more subtle form. We 
show how one can construct an equivalent effective action T[A, D] that is also a functional 
of a classical dimer field D(r,t). This effective action provides a straightforward description 
of states containing a dimer condensate but no atom condensate. 

In section II, we introduce the 1PI effective action for a quantum field theory. We 
determine the 1-body and 2-body terms r 1+2 [A] in the 1PI effective action for the zero- 
range model. In this model, the only interaction is a contact interaction that gives the large 
scattering length a. Any other model with a large scattering length should reduce to the 
zero-range model when \a\ is sufficiently large. In the case a > 0, we construct the equivalent 
terms ri +2 L4, D] in an effective action with a classical dimer field. In section III, we consider 
the approximation defined by the truncated 1PI effective actions r^LA] and r^LA, D\. 
This truncated 1PI approximation corresponds to treating 2-atom quantum effects exactly, 
but ignoring irreducible 3-atom and higher iV-atom effects. To demonstrate the limitations 
of this approximation, we use it to calculate the atom-dimer scattering cross-section. The 
results are compared with those from exact solutions of the 3-body problem, which depend 
on the Efimov parameter k*. In section IV and V, we demonstrate the equivalence of 
the truncated 1PI effective actions by using them to study the strongly-interacting Bose 
gas with total atom number density n at zero temperature. In section IV, we study the 
static homogeneous states. There is an atom superfluid phase and a dimer superfluid phase 
separated by quantum phase transitions at a = and a = (327rn) -1 / 3 . In section V, we 
study the quasiparticles associated with fluctuations around the static homogeneous states. 
The two truncated 1PI effective actions give identical results for the atom superfluid phase, 
but D] is much more convenient for describing the dimer superfluid phase. In section 

VI, we summarize our results and discuss their implications for the Bose gas near a Feshbach 
resonance. 
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II. 1PI EFFECTIVE ACTION 



In this section, we introduce the 1PI effective action for a quantum field theory of bosons. 
If there is a single atom quantum field, the 1PI effective action r[A] is a functional of a 
classical atom field A that encodes complete information on quantum effects. The one-body 
and two-body terms r 1+2 [^4] are given exactly for the zero-range model. For the case a > 0, 
we construct an equivalent functional r 1+2 LA, D] of both the classical atom field A and a 
classical dimer field D. We also show how thermodynamic properties of a many-body system 
can be calculated using these 1PI effective actions. 



A. Quantum field theory 

The simplest quantum mechanical model that can describe atoms with a large scattering 
length a is the zero-range model. In the 2-body sector, the zero-range model can be defined 
by specifying the T-matrix element for atom-atom scattering. If the wavevectors of the two 

— * — * — * — * 

atoms in the initial and final states are ±k and ±k' with \ k\ = \k'\ = k, the T-matrix element 
is 

K J m 1 + iak ' V ; 

where a is the scattering length. Since the T-matrix element is independent of the scattering 
angle, the zero-range model has >S-wave scattering only. The zero-range model plays a 
special role as a minimal model for atoms with large scattering length. In any model of 
atoms whose scattering length a is large compared to the range, the observables at energies 
near the scattering threshold for atoms must approach those of the zero-range model as \a\ 
increases. 

The zero-range model can be formulated as a local quantum field theory. The quantum 
field operator satisfies equal-time commutation relations: 

ft(r, t)ftv\ t) - ftr', t)ft(r, t) = 5(r - r'), (3a) 
ftr,t)ftr',t)-ftv',t)ftr,t) = 0. (3b) 

The Hamiltonian H = J d 3 r7i is the integral of a Hamiltonian density: 

n = ft (-^V 2 ) $ + \\(A)ftftftj>. (4) 
\ 2m ) 4 

The time evolution equation is a partial differential equation: 

+ |^V 2 ) ^(r, t) + l\(A)ft(v, t)ftr, t)ftr, t) = 0. (5) 

The product of operators at the same point in the interaction term is singular, but it can be 
made well-defined by imposing an ultraviolet cutoff that restricts the wavevectors k in the 
Fourier expansion of ip(r,t) to satisfy |k| < A. To have non-trivial scattering in the limit 
A — > oo, the bare coupling constant A(A) must depend on the ultraviolet cutoff A: 

w . . StiU a 

A < A > = —T^r, < 6) 
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Given this coupling constant, the T- matrix element reduces to Eq. (|2j) in the limit A — > oo. 

The model defined by the Hamiltonian density in Eq. is renormalizable in the 2- 
atom sector. This means that all 2-atom observables are well-defined functions of a in 
the limit A — > oo. For example, the model implies the existence of a dimer with binding 
energy Ejj = Ti 2 /{ma 2 ). The model is not renormalizable in the 3-atom sector. To make it 
renormalizable , it is necessary to add a 3-body interaction term (-0 3 )^ 3 to the Hamiltonian 



density |18|, |19J. Its coefficient A3 (A) can be tuned as a function of A so that the Efimov 
parameter has the desired value in the limit A — > 00. All other 3-atom observables are 
well-defined functions of a and in this limit. There is numerical evidence that this 3-b ody 



interaction is also sufficient to make the model renormalizable in the 4-atom sector [20 
If this is true, it is plausible that the zero-range model with a 3-body interaction is also 
renormalizable in the iV-atom sector for all N. 

The model that was taken as the starting point in Refs. 10, III 12|, 13] was a local quan- 
tum field theory with an atom field Vu(r, t) and a molecule field ipM^, t). The Hamiltonian 
density is 

n = ^ A (-|^v 2 ) $ A + A (-|^v 2 + ma)) i> M 

+ \am{^)^m^\^a^m + ^Amm(A) ^ m ^ m ^ m ^m- (7) 

This model is renormalizable in the 2-atom sector where the only parameters are ^m(A), 
A J 4^(A) and #am(A) [2l|. These parameters can be tuned as functions of the ultraviolet 
cutoff A so that the T-matrix element for atom-atom scattering in the limit A — > 00 is 

nk) = -U + - I -f ) \ (8) 

V h k z /m — v J 

where A, g, and v are renormalized parameters that do not depend on the ultraviolet cutoff. 
The scattering length is 

(9) 



m 



8tt?i \ v 




It is not known whether the model defined by the Hamiltonian density in Eq. (JTj) is renor- 
malizable in the 3-atom or higher iV-atom sector. The scattering length in Eq. (Q can be 
made arbitrarily large by tuning to the resonance v — > 0. When this scattering length is suf- 
ficiently large, observables near the scattering threshold for atoms must coincide with those 
of the zero-range model. The scattering length must be large compared to the background 
scattering length mX/(8Trh) and also large compared to the effective range at the resonance, 
which is —lQTch 3 /(m 2 g 2 ). 



B. 1PI effective action 

The one-particle-irreducible (1PI) effective action is a powerful tool for studying quantum 
field theories. The 1PI effective action is introduced in any modern textbook on quantum 
field theory, such as Ref. |22|. In a quantum field theory with only bosonic fields, the 
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1PI effective action is a functional of a set of classical fields, one for each quantum field. 
The classical field associated with a quantum field ip{r,t) that annihilates an atom is a 
complex- valued field that we will denote by A(r,t). We will refer to the field A(r,t) as the 
classical atom field. The 1PI effective action is a functional of A(r,t) that we denote by 
T[A\. Although A(t, t) is a classical field, the functional T[A] encodes complete information 
about the quantum system. 

The 1PI effective action is the generator of 1PI Green functions. The term in T[A] that 
is N th order in both A and A* gives the 1PI iV-body Green function. The 1PI iV-body 
Green functions for n — 1, 2, ■ • • ,N are the building blocks for the iV-body Green function, 
which determines the amplitude for any iV-atom quantum state to evolve into any other 
TV-atom quantum state. The iV-body Green function can be expressed diagrammatically as 
the sum of all possible Feynman diagrams with N lines coming in, N lines going out, and 
arbitrarily many loops. Using the 1PI effective action formalism, the connected iV-body 
Green function can be expressed more compactly as a sum of tree diagrams (i.e. diagrams 
with no loops) whose vertices are 1PI iV-body Green functions with n = 2, 3, . . . , N. The 
1PI Green functions can be obtained from variational derivatives of T[A]. Instead of treating 
the real and imaginary parts of A(r,t) as independent variables, it is more convenient to 
treat A(r,t) and its complex conjugate A*(r,t) as independent. The iV-body 1PI Green 
function can be obtained from T[A] by taking N variational derivatives with respect to A 
and N variational derivatives with respect to A* and then setting A(r, t) = 0. 

The 1PI effective action T[A] also encodes information about the quantum field theory in 
the presence of an external source J(r, t). The zero-range model in the presence of an external 
source is defined by replacing the zero on the right side of the time evolution equation in 
Eq. §5§ by J(r, t). The source is assumed to vanish in the infinite past: J(r,t) — > as 
t — ► — oo. If the system is in the vacuum state in the infinite past, the state will evolve 
under the influence of the external source. In the vacuum state, the expectation value of 
the quantum field vanishes: (ip(r,t)) = 0. In the presence of an external source J(r,t), the 
expectation value of the quantum field can be nonzero and it defines a classical field A(r, t): 

(j>(T 1 t))j = A(T,t). (10) 

The subscript J on the angular brackets in Eq. (flUj) indicates that the expectation value is 
evaluated in the state that evolves from the vacuum under the influence of the external source 
J(r,t). An important basic property of the 1PI effective action is that the first variational 
derivative of T[A] with respect to A*(r,t) gives the source J(r,t) for which Eq. (fTUl) is 
satisfied: 

T[A] = J(r,t). (11) 



6A*(r,t) 



Higher variational derivatives with respect to A(r, t) and A*(r, t) give the 1PI Green function 
for the quantum field theory in the presence of the external source J(r,t). 
If a classical field A(r, t) satisfies the variational equation 

T[A] = 0, (12) 



6A*(T,t) 



Eq. (II ip implies that the corresponding source J(r, t) vanishes. One can then infer from 
Eq. ffTOl) that A(r, t) can be expressed as an expectation value of ip(r,t) in the absence of 
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source. This implies that there is a state \Xa) in the quantum field theory in which the 
expectation value of the quantum field operator is A(r,t): 

(X A ${T,t)\X A )=A{r,t). (13) 

If the source J(r, t) is nonzero but vanishes for all t > t , the vacuum state will have evolved 
at time to m t° a nontrivial state \X A ) of the quantum field theory. The classical field A(r, t) 
given by Eq. (1131) will then satisfy the variational equation (j!2p for all t > t . Thus each 
nontrivial solution A(r, t) of Eq. (1121) corresponds to a nontrivial state \X A ) of the quantum 
field theory. 

The 1PI effective action T[A] can be used to study the static homogeneous states of the 
quantum field theory. If the variational equation ffl2|) has a constant solution A(r,t) = A, 
there must be a static homogeneous state \A) in the quantum field theory in which the 
expectation value of the quantum field is A: 

(A$(r,t)\A) = A. (14) 

The dynamical stability of that static homogeneous state can be studied by considering 
small-amplitude fluctuations of A(r, t) around the constant solution A, i.e. by linearizing 
the variational equation ffl2|) around A. If there are any modes whose amplitudes grow 
exponentially in time, then the state \A) is unstable with respect to quantum fluctuations. 
If there are no such modes, then the state \A) is stable with respect to quantum fluctuations. 

Since it is formulated in terms of the expectation value of a quantum field, the 1PI effective 
action T[A] bears a superficial resemblance to a mean- field approximation. However it is 
not a mean-field approximation, but an exact formulation of the quantum field theory. The 
mean-field approximation for the zero-range model is obtained by replacing the quantum 
field ip{r,t) in the evolution equation (jSJ) by a classical field A(r,t) and by taking the limit 
A — > in the expression for the bare coupling constant A(A) in Eq. (J6]). The resulting partial 
differential equation is 

(. d h 2 „ 9 \ , , . AixTia „, . . . , 

at 2m J 77i 

This equation can be obtained as the variational equation for the classical action S[A\: 

S[ A] = J*j Sr («| + ^V) A - ^(A'Af) . (16) 

Since we have taken the limit A — > 0, the mean-field approximation ignores all effects of 
quantum fluctuations. In contrast, all quantum fluctuations are taken into account exactly 
in the 1PI effective action. 



C. Few-body terms 

If the A^-body Green functions are known analytically, they can be used to deduce the 
A^-body terms in the 1PI effective action. In the zero-range model, the 1-body and 2-body 
Green functions are known analytically, limited analytic information is known about the 
3-body sector 14| , and almost no analytic information is known about the 4-body sector. 
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E — p 2 /2m + ie 
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FIG. 1: Feynman rules for the zero-range model. 
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FIG. 2: The series of Feynman diagrams in the zero-range model whose sum is the 2-atom Green 
function. 



The 1-body term in V[A] can be deduced from the propagator for an atom, which is given 
in Fig. [U The 1-body term is 



r x [A] = J dt J d 3 rA* (ifj + ±V 2 ) A. 



(17) 



Here and in much of the remainder of this paper, we set H = m = 1 for simplicity. Dimen- 
sional analysis can be used to reintroduce factors of h and m when desired. The 1-body term 
in Eq. (|T7|) implies that the energy- momentum relation for an isolated atom is E — p 2 /2m. 

The 2-body term in T[A] can be deduced from the Green function for atom-atom scat- 
tering. The Feynman rules for the atom-atom interaction vertex is shown in Fig. [IJ The 
Green function for atom-atom scattering can be expressed as the infinite sum of Feynman 
diagrams in Fig. [2J Renormalization is accomplished by using Eq. (jBJ) to eliminate A (A) 
in favor of a. The 1PI Green function depends only on the total energy E and the total 
momenta p of the two atoms. Galilean invariance requires it to depend only on E — p 2 /4m. 
The 1PI Green function is 



where the function f(x) is 



A(E,p) = -8vra / (E - p 2 / Am) , 



/(*) 



1 — ay—x — 



IE 



(18) 



(19) 



The T-matrix element for the scattering of two atoms with momenta ±k is obtained by 
evaluating A(E, 0) at the energy E = 2(k 2 /2). After inserting the factors of h and m 
required by dimensional analysis, we recover the T-matrix element in Eq. (J2J). From the 
Green function in Eq. f[T5|) . we can deduce the 2-body term in the 1PI effective action: 



T 2 [A] = -2na J dt J d 3 rA* 2 f(i^ + |V 2 ) A 2 



(20) 



8 



E — p 2 /2m + it 




FIG. 3: Feynman rules for the truncated 1PI effective action ri+2^]- 

The 1PI Green function in Eq. (1181) contains information about bound states in the 
zero-range model. The poles of the amplitude (j!8p on the physical sheet of the complex 
energy E are the energies of two-body bound states. The function f{x) in (|T9"|) has a pole 
at x = — 1/a 2 + ie. If a < 0, there are no bound states, because the pole in E is on the 
unphysical sheet. If a > 0, there is precisely one bound state and its binding energy is 
Ed = 1/a 2 . We refer to this bound state as the dimer. The functional T 2 [A] in Eq. (120]) 
knows about the existence of the dimer through the pole at x = —1/a 2 of the function f(x). 

The 1PI effective action provides a compact representation of iV-body Green functions 
as sums of tree diagrams whose vertices are 1PI Green functions. The Feynman rules for 
the 1-body and 2-body terms in the 1PI effective action are shown in Fig. [31 The Green 
function for atom-atom scattering is given simply by the vertex in Fig. [21 which depends 
only on the total energy E and total momentum p of the two incoming atoms. The T-matrix 
element for the scattering of atoms with momenta ±k is obtained by setting E = 2(k 2 /2m) 
and p = in that vertex. After inserting the factors of % and m required by dimensional 
analysis, this reproduces the T-matrix element in Eq. ([2"1) . 



D. Introducing a classical dimer field 

In the case a > 0, a pair of atoms can bind to form the weakly-bound dimer with binding 
energy Ed = 1/a 2 . The rate for processes involving dimers can be calculated using standard 



diagrammatic methods by using the Lehmann-Symanzik-Zimmermann (LSZ) formalism [22 
This method exploits the fact that the local composite operator ip(r,t) 2 has a nonzero 
amplitude for annihilating a dimer. The rate for processes involving dimers can also be 
calculated using the 1PI effective action by exploiting the fact that 1PI amplitudes have 
poles in the total energy E of any pair of atoms. For example, the atom-atom scattering 
vertex in Fig. [3] has a pole at E — p 2 / (4m) — l/(ma 2 ), which is the energy of a dimer with 
momentum p. One might hope to be able to calculate such processes more conveniently 
using an effective action that is a functional of a classical dimer field D(r,t) as well as the 
classical atom field A(r, t). This motivates the attempt to construct a functional T[A, D] 
that is equivalent to T[A], 

To construct an effective action T[A, D] equivalent to T[A], we require it to generate the 
same variational equation for the classical atom field A after eliminating D. The variational 
equation obtained by varying the functional consisting of the sum of the 1-atom 
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term T^A] in Eq. (HTJ) and the 2-atom term T 2 [A] in Eq. ([20} is 



(i% + ±V 2 ) A = 4vraA7(4 + |V 2 ) A 2 , (21) 

where fix) is the function defined in Eq. (|T9j) . Since /(x) has a pole at x = — 1/a 2 , the 
right side of Eq. f[2~Tj) includes a term in which the operator (i^ + ^V 2 + 1/a 2 ) acts on 
A 2 . We can formally define a classical dimer field D(r,t) by the equation 

D(v, t) = \g (i% + ±V 2 + 1/a 2 ) ^ A 2 (r, t), (22) 

where g is for now an arbitrary constant. The variational equation (J2TJ) can be expressed as 
a pair of coupled equations for A and D: 

(iff + f V 2 ) A = ^^-A*D + W*/(4 + IV 2 ) A 2 , (23a) 

(»* + I V2 ) D = + \9#> (23b) 

where /(#) is a function that is regular at x = —1/a 2 : 

f{x) = f{x) - — 2 2 (24a) 
1 + a l x + xe 



-1 



1 + ay/—x — ie 



(24b) 



The definition of the classical dimer field in Eq. (|22|) follows from Eq. (123bl) . Upon using 
Eq. (1221) to eliminate D from Eq. (I23al) . we recover Eq. (12 ip for the time evolution of the 
classical atom field. 

The Eqs. (1231) can be obtained as variational equations from an effective action T 1+2 [A, D] 
that is a functional of A(r, t) and -D(r, t). The constant g in the definition of D in Eq. (I2"2"j) 
can be determined by demanding that the integrand of the functional includes the standard 
1-body term for a dimer field: D* + |V 2 ) D. The resulting value of g is 

g = (lGTr/a) 1 / 2 . (25) 

The effective action is 

r 1+2 [A, D] = j dt j d 3 r {A* (i| + |V 2 ) A + D* (z| + ^V 2 + 1/a 2 ) D 

-\g iA* 2 D + D*A 2 ) - 2TxaA* 2 fi%% + ^V 2 ) A 2 } . (26) 

By varying this functional, we obtain the pair of coupled variational equations ( 1231) . 

The method used to construct the functional r 1+2 [A,D] in Eq. (|26|) . which includes all 
1-atom and 2-atom terms, could in principle be used to construct a complete 1PI effective 
action T[A, D] that includes all higher iV-atom terms. If the iV-atom Green function were 
known analytically, terms that contain poles in the total energy of an incoming or outgoing 
pair of atoms would be replaced by terms with a classical dimer field D or D*. The resulting 
1PI effective action T[A, D] for classical atom and dimer fields would be equivalent to the 
conventional 1PI effective action T[A\. 
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E — p 2 /2m + it 

E — p 2 /4m + 1/ma 2 + it 




FIG. 4: Feynman rules for the truncated 1PI effective action ri + 2[A-D]- 

Using the 1PI effective action T[A, D\, the Green function for N atoms could be repre- 
sented as the sum of tree diagrams with N incoming atom lines, N outgoing atom lines, and 
vertices connected by either atom lines or dimer lines. The Feynman rules for the 1-atom 
and 2-atom terms in the effective action r 1+2 [A, D] are shown in Fig. HI They include prop- 
agators for the atom and the dimer, an atom-atom scattering vertex, and a point coupling 
between two atoms and a dimer. The Feynman rule for the atom-atom scattering vertex in 
Fig. H] is different from the one in Fig. [3J The Green function for atom-atom scattering is 
the sum of the two tree diagrams in Fig. The T-matrix for the scattering of atoms with 
momenta ±fc is obtained by setting E = 2(k 2 /2m) and p — 0. After inserting the factor of 
h required by dimensional analysis, this reproduces the T-matrix element in Eq. (J2]). 

The 1PI effective action r[A, D] also provides a convenient representation of amplitudes 
for processes with dimers in the initial and final states. Such an amplitude is given by the 
sum of tree diagrams with an external dimer line for each of the initial or final dimers. This 
amplitude could also be derived from T[A] by using the LSZ formalism [22| which requires 
finding the residues of poles in the total energy of pairs of external atom lines. The effective 
action T[A, D] gives these amplitudes much more directly. 

The value of the atom-dimer coupling constant g in Eq. (1251) has been obtained previously 
in Ref. The authors expressed the atom-molecule interaction term in the Hamiltonian 
density Eq. as \/~Zg { r tp' [ M 'ipA'^ l A+h.c.) , where \/~Z is the amplitude for a "dressed molecule" 
to be a "bare molecule". As v — > 0, which implies a — > ±oo, the limiting behavior of this 
amplitude is \J~Z — > 2y/nh 2 /(gmy^a). Setting % — m — 1, the limiting behavior of the 
coefficient of ^ M ipAipA is V~Zg — > a/47t /a. This agrees with the coefficient \g of the term 
D*A 2 in the integrand of the effective action functional in Eq. fl26]) . 
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FIG. 5: Feynman diagrams for atom-atom scattering from the truncated 1PI effective action 
T[A,D}. 

E. Many-body problem 

The 1PI effective action can also be applied to a many-body system containing infinitely 
many atoms with a given local number density n(r, t). IfH is the Hamiltonian density oper- 
ator, the many-body system can be described by a quantum field theory whose Hamiltonian 
density operator is Ti. — \i ip^ip, where \i is the chemical potential. In the iV-body sector, the 
only effect of a constant chemical potential \x is to shift the total energy by Nfj,. If the system 
is homogeneous, a can be adjusted so that the number density n has the desired value. In a 
nonhomogeneous many-body system, the function /i(r,t) can be chosen so that the ground 
state of the system has the desired local number density n(r, t). Position-dependent terms 
in /i(r, t) play the rule of a trapping potential. In the case of a time- independent chemical 
potential /x(r), the 1PI effective action r[A;/x] can be obtained from the conventional 1PI 
effective action T[A] simply by replacing the time derivative i-j^ with i-jL + /j when it acts on 
A as in Eq. (1T7I) . with i-^ + 2n when it acts on A 2 as in Eq. fl2Dj), etc. If the chemical potential 
jti(r) is time dependent, obtaining T[A; /i] is more difficult because \i does not commute with 
the time derivative J^. 

The local number density n(r, t) is the expectation value of the number density operator: 

n(r,t) = (^(r,t)^(r,t)>. (27) 

It is the conserved density associated with the invariance of the 1PI effective action T[A; fj] 
under the phase symmetry A — > e ia A. It can therefore be derived from the effective action 
using the generalization of Noether's prescription to the case of arbitrarily high-order time 
derivatives. Alternatively, it can be obtained simply by varying T[A; jj\ with respect to the 
chemical potential: 

n ^ ()= h^) miA (28) 

The 1PI effective action r[A, D] with a classical dimer field can also be generalized to a 
functional T[A, D; [i] that describes the many-body system with chemical potential n(r,t). 
From the definition of the classical dimer field in Eq. (1221) , we can see that its behavior under 
the phase symmetry is D — > e t2a D. In the case of a time-independent chemical potential 
/x(r), T[A,D;/j] can be obtained from T[A, D] by replacing the time derivative ij^ with 
zj^ + fx when it acts on A, with + 2ji when it acts on A 2 or D, etc. The number density 
n(r, t) can be determined by using the fact that it is the conserved density associated with 
the phase symmetry. Thus it can be obtained by varying T[A, D; fi] with respect to the 
chemical potential as in Eq. fl28|) . 

The energy density £ is the expectation value of the hamiltonian density operator: £ = 
(H). The free energy density T = £ — fin is the conserved density associated with the 
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invariance of the effective action T[A; /i] under translations in time: A(r, t) — > A(r, t + e). It 
can therefore be derived by using Noether's method generalized to the case of arbitrarily high 
order time derivatives. Similarly, if we use the effective action T[A, D; /j] with a classical 
dimer field, the free energy density can be determined by using the fact that it is the 
conserved density associated with translations in time of A(r,t) and D(r,t). 



III. TRUNCATED 1PI EFFECTIVE ACTION 

In this section, we introduce an approximation in which the 1-body and 2-body terms 
in the 1PI effective action are included exactly, but 3-body and higher-body terms are 
neglected. We refer to this approximation as the truncated 1PI approximation. The cross 
section for atom-dimer scattering is calculated using the truncated 1PI approximation and 
compared to results from the exact solution of the 3-body problem. 



A. Classical atom field only 

We define the truncated 1PI effective action for the many-body system by 

T 1+2 [A;fi] = J dt J d\{A* (t^ t +ii+\V 2 )A-27raA* 2 f(t^ t +2fi+\V 2 )A 2 }, (29) 

where f(x) is the function defined in Eq. (fl9l) . The variational equation for the classical 
field A(r,t) is 

{%% + fi + ±V 2 ) A = 47raA7(z| + 2^ + ±V 2 ) A 2 . (30) 

Note that this equation is nonlocal in position and time, because / is a nonpolynomial 
function of an operator that involves the derivatives and V. A solution of this equa- 
tion corresponds to a state containing a Bose-Einstein condensate of atoms whose common 
wavefunction is proportional to A(r, t). 

The number density n(r, t) that follows from the truncated 1PI effective action in Eq. (f29|) 
is given by the variational derivative with respect to fi(r,t) in Eq. (f28|) . In the special case 
of a constant chemical potential /i, the variational derivative can be evaluated explicitly: 

n = A* A - AnaA* 2 f (i^ + 2/i + ±V 2 ) A 2 , (31) 

where f'(x) is the derivative of the function defined in Eq. (|T9|) : 

/'(*) = ~ Q , T2- (32) 

2v— x — it (1 — ay—x — le) 

The free energy density !F(r, t) that follows from the truncated 1PI effective action in 
Eq. (1291) can be obtained by applying Noether's theorem to the time-translation symmetry. 
In the special case of a time- independent classical field A(r), Noether's prescription reduces 
to 

T = -A* (fi + \ V 2 ) A + 27raA* 2 f (2/i + \ V 2 ) A 2 . (33) 
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B. Classical atom and dimer fields 



In the case a > 0, we can construct an equivalent truncated 1PI effective action that is a 
functional of classical atom and dimer fields: 

T 1+2 [A, D;n] = J dt J d 3 r {A* (ifj + /i + |V 2 ) A + D* (i% + 2/i + |V 2 + 1/a 2 ) D 

-iff (A* 2 D + D*A 2 ) - 27caA* 2 f(t^- t + 2// + |V 2 ) A 2 } , (34) 

where / (x) is defined in Eq. (I24p and g = ^lQir/a. The equations obtained by varying this 
functional with respect to A*(r,t) and D*(r,t) are 

(«ft + M + |V 2 ) A = — A*D + 47raA7(i|j + 2/i + ±V 2 ) A 2 , (35a) 

aff 

(i| + 2/i + iV 2 ) -D = -Id + \gA 2 . (35b) 
The formal solution of Eq. (135bj) for the field D(r,t) is 

D(r, t) = \g (i| + 2^ + J V 2 + 1/a 2 ) _1 A 2 (r, t). (36) 

Upon using this solution to eliminate D from Eqs. (I35ap . we recover Eq. (|30|) for the classical 
atom field. A solution of Eqs. fl35l) corresponds to a state containing a condensate of atoms 
whose common wavefunction is proportional to A(r, t) and a condensate of dimers whose 
common center-of-mass wavefunction is proportional to D(r,t). 

The number density n(r,t) associated with the truncated 1PI effective action in Eq. (13"4"|) 
is given by the variational derivative with respect to fi(r, t) in Eq. f[2"5|) . In the special case 
of a constant chemical potential /i, the variational derivative can be evaluated explicitly: 

n = A* A + 2D*D - AnaA* 2 /' (i^ + 2/i + \ V 2 ) A 2 , (37) 
where f'{x) is the derivative of the function defined in Eq. (J21 



m = — - a/2 . 2/ . . (ss) 

(1 + a\J —x — 2e) z V — ^ — *e 



The energy density £(r, t) that follows from the truncated 1PI effective action in Eq. (jH! 
can be obtained by applying Noether's theorem to the time-translation symmetry. In the 
special case of time-independent classical fields A(r) and D(r), Noether's prescription re- 
duces to 

T = -A* (fi + f V 2 ) A — D* (2/i + iV 2 + 1/a 2 ) D 

+ \g (A* 2 D + D*A 2 ) + 2vraA* 2 /(2/i + |V 2 ) A 2 . (39) 



C. Atom-dimer elastic scattering 

The truncation of the 1PI effective action to include only 1-atom and 2-atom terms 
guarantees that 2-atom effects are treated exactly. In this approximation, irreducible N- 
atom amplitude for N > 3 are neglected. This does not imply that higher A-atom effects 
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FIG. 6: Feynman diagram for elastic atom-dimer scattering from the truncated 1PI effective action 

r 1+3 [A,z>]. 



are trivial, because they can arise from the iteration of 2-atom effects. In this subsection, 
we clarify the meaning of the truncated 1PI effective action by using it to calculate the cross 
section for atom-dimer elastic scattering. We will compare the cross section with results 
from exact solutions of the 3-body problem. In this subsection, we set \x = and make the 
mass m explicit in the equations instead of setting m = 1. 

If a > 0, the truncated 1PI effective action r 1+2 [A] gives nontrivial predictions for atom- 
dimer scattering that can be obtained by using the LSZ formalism. These predictions can be 
obtained more directly from the equivalent truncated effective action Ti + 2[A D] in Eq. fl34|) 
with a classical dimer field. Atom-dimer scattering proceeds through the Feynman diagram 
in Fig. [61 The T-matrix element obtained using the Feynman rules in Fig. H] is 

T(k 9) « 16na/m (40) 

i(/e,tf)~ i + ( 5+4cos 0) a2A .2/ 4 > l 4U J 

where k is the momentum of the atom or the dimer in the center-of-mass frame and 9 is the 
scattering angle. The atom-dimer scattering length a ad can be obtained from the T-matrix 
element at k — 0. Using T(k, 9) — > —'ina^jm as k — > 0, we obtain 

a-AD = — g-a- (41) 



The atom-dimer scattering length has been calculated previously in Ref. [llj using an ap- 
proximation equivalent to our truncated 1PI approximation. Their result 32a/3 differs from 
ours by a factor of —2. Part of the discrepancy seems to come from a factor of 2 error in 
the expression for the T-matrix element T am associated with the diagram in Fig. 2(a) of 
Ref. [ll| , which is the same as our Fig. [6j 

Using the T-matrix element in Eq. (HUj) . we can obtain the total elastic atom-dimer cross 
section: 

1 024vra 2 /9 

aM ' ~ (l + 3ma 2 £)(l + ma 2 T/3)' { ' 

where E = 3k 2 /Am is the collision energy. By applying the partial wave decomposition 
to the T-matrix element in Eq. (1401) . we can isolate the S"-wave contribution to this cross 
section: 

(£=%e* 1R 2 r 1 , ( 1 + 3ma 2 T 



_ma 2 E ^ \1 + ma 2 E/3j\ ' ^ 
At the dimer-breakup threshold E = Ed = l/(ma 2 ), the S- wave cross section is 

cri L D =0) (T D ) = (16 log 2 3) vra 2 w 19.31 vra 2 . (44) 
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The result in Eq. (I42p for the total cross section at this energy is (64/3)7ra 2 . Thus the 
truncated 1PI effective action predicts that S-wave scattering should account for 90.5% of 
the cross section at the dimer-breakup threshold. 

The results from the truncated 1PI effective action can be compared with those from 
exact solutions of the 3-body problem. The exact results that are known are summarized in 



Ref. [14j. The exact result for the atom-dimer scattering length is 

a ad = (1-46 - 2.15 tan [s ln(aK*) + 1.16] ) a, (45) 

where So = 1.00624 and is the Efimov parameter. The functional form was first de- 
duced by Efimov [l7j]. The numerical coefficients were first calculated by Simenog and 



Sitnichenko [23j and by Bedaque, Hammer, and van Kolck [18|, [l9j. The ratio ciad/o from 
Eq. f )45l) can have any value between — oo and +oo depending on the value of the product 
an*. Thus the result for a ad in Eq. (j4Tl) from the truncated 1PI effective action can have 
an arbitrarily large error. 

Another exact result for atom-dimer scattering is the S-wave phase shift at the dimer- 



breakup threshold, which was derived by Macek, Ovchinnikov, and Gasaneo [24|. The 
S'-wave contribution to the cross section at the dimer-breakup threshold is 

af D =0) (£ D ) = 3rta 2 sin 2 [s ln(a^) + 2.91] . (46) 

The exact cross section can have any value between and 3ira 2 depending on the value of the 
product aft*. The S'-wave cross section from the truncated 1PI effective action in Eq. (14"4")) 
is about 6 times larger than the maximum value of the exact cross section. We see again 
that the truncated 1PI effective action can give results for low-energy 3-body processes with 
large errors. 



IV. STATIC HOMOGENEOUS STATES 

In this section and the next one, we use the truncated 1PI approximation as a model for 
a Bose gas of atoms with large scattering length a and total atom number density n at zero 
temperature. For a > 0, the results illustrate the equivalence of the truncated 1PI effective 
actions Ti + 2[^4; /i] and r^LA, D; /j]. In this section, we study the static homogeneous states 
predicted by this model. 



A. Atom condensate 

We can study the static homogeneous states of the system using the truncated 1PI effec- 
tive action Pi + 2[^4; fi] in Eq. (1291) . In a static homogeneous state, the classical atom field is 
a complex constant: A(r, t) = A. We will refer to A as the atom mean field. A state with 
A ^ contains a Bose-Einstein condensate of atoms. In such a state, a macroscopic fraction 
\A\ 2 /n of the atoms have the same constant wavefunction whose phase is that of A. 

The atom mean field A must satisfy Eq. (13"U|) . which reduces to 

fI A = —\A\ 2 A. (47) 

1 — r 
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We have introduced a dimensionless chemical potential variable r defined by 



r = a(-2fi - ie) 1/2 . (48) 

We assume that \x is negative so r is real-valued and has the same sign as a. The expression 
for the chemical potential in terms of r is 

fi = -r 2 /(2a 2 ). (49) 

If Eq. (T47j) has a nonzero solution A, the atom mean field must satisfy 

\ A \ =— 5-^- 50 

Since the left side of Eq. f )50|) is real and positive, the right side must also be real and 
positive. This requires r < 1 if a < and r > 1 if a > 0. However, if a < 0, Eq. fT4H|) with 
H < sets a stronger constraint r < 0. Thus a static homogeneous state with non-zero A is 
possible only if r < or r > 1. 

The number density n and the free energy density T in the static homogeneous state are 
given by setting A(r,t) = A in Eqs. (jHTl) and (1551) : 

r[r — l) z 

T = ^\A?-^\A\\ (51b) 
2a z r — 1 

Using Eq. (1501) to eliminate \A\ 2 , we obtain parametric expressions for n and as functions 
of r: 

r 2 (5r — 4) 

n = \ 3 J , 52a 

If a < 0, the number density n increases from to oo as r decreases from to — oo. If 
a > 0, the condition r > 1 implies that n is greater than the critical number density 

n c = l/(327ra 3 ). (53) 

Equivalently it implies that a is greater than the critical scattering length 

a c = (32nn)- 1/3 . (54) 

The number density in Eq. (I52al) increases from n c to oo as r increases from 1 to oo. 

The two terms in n and in T in Eqs. ( 1511) can be interpreted as condensate and non- 
condensate contributions, respectively. Dividing Eq. (1501) by Eq. (I52aj) . we obtain the fraction 
of the atoms that are in the atom condensate: 
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Atom + Dimer 




FIG. 7: Energy per particle £/n in units of 1/a 2 as functions of kpa, where kp = (67r 2 n) 1//3 . The 
vertical dotted lines separate the regions of kpa in which the static homogeneous state includes 
an atom condensate (a < 0), a dimer condensate (0 < a < a c ), and a mixture of atom and dimer 
condensates (a > a c ). The dashed line represents an unstable homogeneous state that includes an 
atom condensate. 



If a < 0, this fraction increases from to 4/5 as n increases from to oo. If a > 0, this 
fraction increases from to 4/5 as n increases from n c to oo. 

The energy density in the static homogeneous state is £ = T + fin, where /x, n, and T 
are given in Eqs. fl4T?|) and (1521) . The average energy per atom is 

e t n = - r ;< 3r - 2 > (56) 
7 2(5r-4)a 2 V ; 

The energy per atom is shown as a function of a in Fig. [3 In the high-density limit, a 
convenient energy scale is the Fermi energy ep for an ideal gas of fermions with a single spin 
component: ep = |(67r 2 n) 2 / 3 . The energy per atom in Eqs. (|56l) has a well-defined unitary 
limit as a —>■ ±oo with n fixed: 

3 / 16 \ 2 / 3 



£/n — j e^coo. (57) 

Its numerical value is eoo ~ —0.292 €f. If a < 0, S/n decreases from 27ma to e ro as n increases 
from to oo. The limiting value 2ima as n —>■ is the analytic continuation to negative a 
of the energy per particle of a weakly interacting Bose gas. If a > 0, £/n decreases from 
— l/(2a 2 ) to as n increases from n c to oo. The limiting value — l/(2a 2 ) at n = n c is the 
energy per atom for a pair of atoms bound into a dimer with binding energy Ep, = 1/a 2 . 
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If a > and < n < n c , there are unphysical solutions to Eqs. fH7|) and fl51al) . The 
expression for n in Eq. (I52ap is a cubic polynomial in r. It therefore has three roots: a real 
root in the range | < r < 1 and a pair of complex roots. For the real root, the expression 
for | A | 2 in Eq. (I50p has an unphysical negative value. For the complex roots, the expression 
for \A\ 2 in Eq. (1501) has an unphysical imaginary part. If n <C n c , the complex roots can be 
expanded in powers of na 3 : 

r = ±i(8nna 3 ) 1/2 [l ± i§(8vma 3 ) 1/2 + ...]. (58) 

The thermodynamic quantities corresponding to this root can also be expanded in powers 
of na 3 . The expansion for the energy per atom is 

£/n = 2nna [l ± z(8vma 3 ) 1/2 + ...]. (59) 

The leading term in Eq. (1591) is the well-known result for the weakly-interacting Bose gas. 
If na 3 is small, the imaginary part is suppressed by a factor of (na 3 ) 1//2 . If a is decreased 
through to a negative value, one of the roots in Eq. (I5"8"j) becomes the negative value of 
r that corresponds to the static homogeneous state with mean- field A satisfying Eq. fl50l) . 
Thus the weakly-interacting Bose gas can be interpreted as the analytic continuation to 
a > of the static homogeneous state for a < 0. The real part of ^/n as a function of a is 
shown clS db dashed line in Fig. [7J 



B. Mixture of atom and dimer condensates 

If a > 0, we can also study the static homogeneous states of the system using the 
truncated 1PI effective action Fi +2 [A, D; /j] given in Eq. flMl) . In a static homogeneous state, 
the classical atom and dimer fields must be complex constants: A(r,t) = A, D(r,t) = D. 
We will refer to A and D as the atom and dimer mean fields, respectively. A state with 
and D ^ contains a mixture of Bose-Einstein condensates of atoms and dimers. In such 
a state, a macroscopic fraction |A| 2 /n of the atoms have the same constant wavefunction 
whose phase is that of A. A macroscopic fraction 2\D\ 2 /n of the atoms are bound into 
dimers that have the same constant center-of-mass wavefunction whose phase is that of D. 

The atom and dimer mean fields A and D must satisfy Eqs. (1351) . which reduce to 

fjJi = —A*D - —\A\ 2 A, (60a) 
ag 1 + r 

2/iD = -^D + \gA 2 , (60b) 

where r is the dimensionless chemical potential variable defined in Eq. ( )48|) . Eq. (160bl) can 
be solved for the mean field D function of A: 

D = - 9 f A 2 . (61) 
2(r 2 - 1) V ; 

Inserting this expression for D into Eq. (I60al) and solving for \A\ 2 , we obtain the same 
expression as in Eq. (1501) . which was derived from the action Fi +2 [v4; //]. 
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The number density n and the free energy density T in the static homogeneous state are 
given by Eqs. (1371) and ([39]) : 

n = \A\ 2 + 2\D\ 2 + 27Ta .j A\\ (62a) 
r(r + iy 

T = ^\ 2 + r ^\Df + 9 Re(D*A 2 )-^\A^ (62b) 

The first two terms in n and in T can be interpreted as atom condensate and dimer 
condensate contributions, respectively. The remaining terms can be interpreted as non- 
condensate contributions. Using Eqs. (ISTj) and (1501) to eliminate the mean fields D and A 
in favor of r and a, we obtain the same expressions for n and T as in Eqs. ( |52|) . which were 
derived from the action ji}. 

Nonzero values of A and D that satisfy the variational equations for Ti + 2[A, D; fj] corre- 
spond to a state that contains a mixture of an atom condensate and a dimer condensate. One 
might be tempted to interpret a nonzero value of A that satisfies the variational equations 
for r 1+2 [A; /j] as the atom mean field in a state with an atom condensate but no dimer con- 
densate. In this case it would seem quite remarkable that the static homogeneous state with 
nonzero A and D derived from D; fj] has exactly the same thermodynamic properties 

as the static homogeneous state with the same value of A derived from ri + 2[A; fi]. However 
the functionals T 1+2 [A; fi] and r 1+2 [A, D; /i] are equivalent. Thus the state with nonzero A 
and D derived from r 1+2 [A, D; yu] is in fact identical to the state with the same value of A 
derived from r 1+2 [A;/i]. Constructing the equivalent action T 1+2 [A, D; fi] and looking for 
solutions with a nonzero value of D as well as a nonzero value of A does not reveal any 
new states of the system. However it does allow non-atom-condensate contributions to ther- 
modynamic quantities to be resolved into dimer condensate contributions and contributions 
that correspond to neither atom nor dimer condensates. 

Using Eqs. (ISTj) . (1501) . and (I52al) . we can determine the fraction of the atoms that are in 
the dimer condensate: 

2\D\ 2 /n=- ^ -. (63) 

This fraction decreases from 1 at n = n r to as n — > oo. 



C. Dimer condensate 

If a > 0, we have found that for number density n less than the critical value n c defined 
in Eq. (1531) . there are no static homogeneous states with A ^ 0. However, there is another 
possibility. If the chemical potential has the value 

/i = -l/(2a 2 ), (64) 

the variational Eqs. (1601) can be satisfied if A — and 5^0. This corresponds to a state 
that has a dimer condensate but no atom condensate. Since r = 1, the number density and 
the free energy density in Eqs. (1621) reduce to 

n = 2\D\ 2 , (65a) 
T = 0. (65b) 
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The energy density is therefore S = jj,n } and the average energy per atom is 

S/n = -l/(2a 2 ). (66) 

The energy per atom in this state is shown as a horizontal line in Fig. [7J Since the number 
density of dimers is n/2, the energy per dimer is — 1/a 2 . Both a and S/n have the same 
values as a gas consisting of dimers that have the binding energy Ed = 1/a 2 but are otherwise 
noninteracting. Thus we can interpret this state as a pure dimer condensate. 



D. Discussion of phases 



The truncated 1PI approximation predicts one static homogeneous state for any value of 
the scattering length a. We can interpret the various static homogeneous states as different 
superfluid phases separated by phase transitions at a = and at the critical value a = a c 
defined in Eq. 

• an atom superfluid phase with only an atom condensate for — oo < a < 0, 

• a dimer superfluid phase with only a dimer condensate for < a < a c , 

• an atom superfluid phase with a mixture of atom and dimer condensates for a c < a < 
+oo. 

The phase transitions are illustrated in Fig. [7J which shows the average energy per atom 
S/n as a function of the dimensionless variable kpa, where kp = (67r 2 n) 1 / 3 is the Fermi 
wavenumber of a system of fermions with a single spin state and number density n. The 
energy per particle is given by Eq. floT)]) for — oo < a < 0, by Eq. f[6"6"l) for < a < a c , and by 
Eq. ([SSD again for a c < a < +oo. The energy per particle is discontinuous at a = and its 
second derivative is discontinuous at a = a c . 

We proceed to determine the prediction of the truncated 1PI effective action for the 
orders of the two phase transitions. The free energy density T is given by Eq. (152bl) for 
— oo < a < 0, by Eq. (165bl) for < a < a c , and by Eq. (152b|) again for a c < a < +oo. The 
limiting behavior of T in the atom superfluid phase as a approaches the phase transition at 
a = is 

T — > —luan 2 as a — > CP. (67) 

Since T = in the dimer superfluid phase, T is continuous at a = but its first derivative 
with respect to a is discontinuous. Thus the phase transition at a = is predicted to be 
first order. The limiting behavior of T in the atom superfluid phase as a approaches the 
phase transition at a = a c is 

T — ► n 2 (a — a c ) as a — >a^. (68) 

Since T = in the dimer superfluid phase, T is continuous at a = a c but its first derivative 
with respect to a is discontinuous. Thus the phase transition at a = a c is also predicted to 
be first order. 

The free energy density in Eq. (I52bl) has the same limiting behavior as a —>■ +oo and 
a — > — oo, so it has a well-behaved unitary limit: 

T — (327r ( 5)2/3 n5/3 as a - ±oo. (69) 
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FIG. 8: Energy per particle £ /n in units of the Fermi energy ef = k F /2 as functions of (/cfci) -1 , 
where kp = (67r 2 n) 1//3 . The vertical dotted lines separate the regions of (fc^a) -1 in which the 
static homogeneous state includes an atom condensate (a -1 < 0), a mixture of atom and dimer 
condensates (0 < a -1 < a^ 1 ), and a dimer condensate (a -1 > a^ 1 ). 

Not only is T a continuous function of a -1 at a -1 = 0, but all derivatives of T with respect 
to a~ l are also continuous. Thus there is no phase transition at a~ l = between the atom 
superfluid phase for a -1 < and the atom superfluid phase for a" 1 > 0. The system in these 
two regions differs in that there is a dimer condensate for a -1 > and no dimer condensate 
for a -1 < 0. However the fraction of atoms in the dimer condensate goes to in the unitary 
limit in such a way that all thermodynamic variables remain smooth. The smooth behavior 
of the energy per atom as a function of a -1 at a^ 1 = is illustrated in Fig. [SJ 

The smooth behavior of the thermodynamic variables at a -1 = is consistent with 
the universality hypothesis proposed by Ho [25], which asserts that the thermodynamic 
properties of the system in the unitary limit a = ±oo are universal functions of the number 
density n and the temperature T and don't depend on any interaction parameters. This 
hypothesis requires that a thermodynamic variable at T = and a = ±oo have the form 
A ri p , where A is a constant and p is the power required by dimensional analysis. This 
hypothesis is well established for systems consisting of fermions with two spin states. Ho 
proposed that it might also apply to bosonic systems. 

If a many-body system consisting of identical bosons with large scattering length a has 
number density n and temperature T such that ri _1//3 and (mT/h 2 )' 1 ^ 2 are large compared to 
the range, its behavior should be determined by a and the Efimov parameter and it should 
be governed by the discrete scaling symmetry. Efimov physics requires a modification of Ho's 
universality hypothesis [25] to allow for log-periodic dependence on the Efimov parameter k*. 
For example, a thermodynamic variable at zero temperature in the unitary limit a = ±oo 
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must have the form n p f(so log(n//^)), where p is the power required by dimensional analysis 
and f{6) is a dimensionless function with period 2ir. Our truncated 1PI approximation does 
not give this behavior, because irreducible 3-body amplitudes have been neglected. 



E. Comparison with previous work 

Several groups have studied the static homogeneous states for the many-body system of 
bosonic atoms near a Feshbach resonance 10, 11, Qjl, HI]. The starting point for all these 



analyses could be taken to be the quantum field theory with an atom field Tp^(r,t) and a 
molecule field if}M{r,t) defined by the Hamiltonian density in Eq. ([7]). They considered the 
phase diagram of the system as a function of the atom number density n, the temperature 
T, and a detuning variable that determines the scattering length a of the atoms. When the 
scattering length is large compared to all other length scales, the low energy behavior of the 
system should be the same as in the zero-range model. We can compare our results from 
the truncated 1PI effective action at zero temperature with their results for large scattering 
length. 

We first consider the qualitative behavior of the phase diagram at T = 0. Radzihovsky, 
Park and Weichman (RPW) and Romans, Duine, Sachdev and Stoof (RDSS) independently 
pointed out that this system has a quantum phase transition between an atom superfluid 



phase and a molecular superfluid phase [10|, 111)]. The quantum phase transition is between 



states with different symmetries. The fundamental theory defined by the Hamiltonian den- 
sity in Eq. ([7]) has a U(l) symmetry corresponding to phase transformations of the atom and 
molecule quantum fields: ipA e^i^A, i>M ~^ e VW- In the atom superfluid phase, both 
the atom and molecule fields have nonzero expectation values (ip a) and (iPm), so the U(l) 
symmetry is completely broken to the trivial subgroup {+1}. In the molecular superfluid 
phase, (V>m) is nonzero but (tpA) = 0, so the U(l) symmetry is broken to the subgroup 
Z 2 = {+1,-1}. Thus the order parameter for the phase transition is an element of the 
discrete group Z 2 [l(J [Hj]. For a > 0, the classical atom and dimer fields A and D in the 
truncated 1PI effective action r 1+2 [A, D; /i] play the same role as the expectation values 
(ipA) and (V>m) in Refs. [l(J U|- For a < 0, we have no analog of (^m)- However the 



classical atom field A in the truncated 1PI effective action Pi +2 [A; /i] plays the same role as 



the expectation value (i/^a) in Refs. [10j, 111 



The order of the quantum phase transition between the atom superfluid phase and the 



molecular superfluid phase has been discussed in Refs. [10|, [III, ll2| . Because the two phases 
are distinguished by a Z 2 order parameter, one might expect the phase transition to be in 
the universality class of the 4-dimensional Ising model, which has a second order transition. 
However, the coupling of the Z 2 order parameter to the Goldstone mode associated with 
the molecule condensate phase modifies the critical behavior and makes the phase transition 
weakly first order. The mean-field approximation predicts that the phase transition is second 
order. Lee and Lee used the renormalization group to show that the phase transition is first 
order The truncated 1PI approximation correctly predicts that the phase transition is 
first order. 

The analysis of the static homogeneous states by RPW [lOj was carried out using the 
mean-field approximation, in which the quantum field operators in Eq. (j7|) are replaced by 
classical fields ^a( t , t) and iPm{y, t) and the bare parameters, which depend on the ultraviolet 
cutoff A, are replaced by mean-field parameters. The analysis of the static homogeneous 
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states by RDSS [11] included some effects that go beyond the mean-field approximation. 
The grand Hamiltonian density in the mean-field approximation is 



H - \m = ij)* A (-±V 2 - //) + ip*M (-|V 2 - 2/x + i/m) 

+ ^AM^*M^P*A^A^M + \X M M {^Ii^mT • (70) 

The mean-field predictions for the scattering observables are sums of tree diagrams whose 
vertices are given by the interaction terms in Eq. (1701) . For example, the T-matrix element 
for atom-atom scattering is given by the sum of the two diagrams in Fig. where the 
vertices in the first diagram are —igAM and the vertex in the second diagram is — i\aa- The 
resulting expression for the scattering length is 

a=^(\AA- 9 ^). (71) 



8tt V. vm 

The resonance where the scattering length a diverges is reached by tuning the parameter 
vm to 0. 

We first discuss the values of the mean-field parameters in Ref. [Io|. RPW stated that 
in the dilute gas limit, Xaa, Xam, Xmm are proportional to the scattering lengths for 
atom-atom, atom-molecule, and molecule-molecule scattering, respectively. This presum- 
ably means that the scattering length a is given by Eq. (17TT) with qam = 0. This would 
be inconsistent with the classical approximation. RPW did not specify the parameter gAM- 
They also did not discuss the behavior of the mean-field parameters near the resonance 
where the scattering length a is large. 

The approximations of RDSS can also be expressed in terms of an effective Hamiltonian 
with the terms in Eq. (I70p [ill ]. They discussed the behavior of the mean-field parameters 
near the resonance where a — > +oo. For a > 0, which corresponds to negative detuning, the 
limiting behavior of the parameters is 

(72a) 
(72b) 
(72c) 
(72d) 
(72e) 

The limiting behavior of the detuning parameter um is the energy of the weakly-bound 
dimer. The limiting value of gAM also comes from the exact solution of the 2-body problem. 
Note that all the coefficients in Eqs. (f72|) scale as powers of a except Xaa- RDSS took 
Xaa to be a constant, so it gives a non-resonant contribution db g to the scattering length 
in Eq. (17TT) . Note that if the limiting values of u M , Xaa, and g A M i n Eqs. (172]) are inserted 
into the right side of Eq. (I7T1) . the limiting value is 2a instead of a. RDSS took Xam and 
Xmm to be proportional to the scattering lengths for atom- molecule and molecule- molecule 
scattering. They recognized that in the limit a — ► oo, these scattering lengths must scale 
linearly with the large scattering length a. RDSS did not recognize that Efimov physics 
requires the coefficients of a in these scattering lengths to be log-periodic functions of a 
of the form /(s log(a/«*)), where is the Efimov parameter and f{6) is a function with 



vm — * 


-1/a 2 , 


gAM — * 




Xaa = 


8ira bg , 


Xam — > 


327T<2, 


Xmm — > 


lQna. 
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period 2ir. They calculated the atom-dimer and dimer-dimer scattering lengths using an 
approximation in which 2-body quantum effects are treated exactly, the atom-dimer coupling 
constant qam is treated as a perturbation, and other irreducible 3-body and 4-body quantum 
effects are neglected. In the case of the atom-dimer scattering length a ad, this approximation 
is equivalent to our truncated 1PI approximation. Their result for a ad differs from ours in 
Eq. (14 ip by a factor of —2. Both results for a ad differ from the exact result in Eq. (1451) in 
which the coefficient of a is a log-periodic function of a/^. The approximation of RDSS for 
the dimer-dimer scattering length is a^D = 4a. The exact result for aoD has not yet been 
calculated. 

Our truncated 1PI approximation gives predictions for the parameters in the mean-field 
Hamiltonian in Eq. (ITUl) in the limit of large scattering length. The classical atom field A(r, t) 
in the 1PI effective action can be naturally identified with the atom mean field iPa{ v , t). If 
a > 0, we can also identify the classical dimer field D(r,t) with the molecular mean field 
ijjM( r ,t)- Inserting constant mean fields A and D into the truncated 1PI effective action in 
Eq. fl3"4"|) . it reduces to 

r 1+2 [A,D;v) = J dt J d 3 r {/j,A*A + (2/i + 1/a 2 ) D*D 

-^fa{A* 2 D + D*A 2 ) + 2na r _ {A* A) 2 ) . (73) 

1 + fly — 2/i J 

We can identify the integrand with the negative of the Hamiltonian density in Eq. ( JTOl) for 
constant values of the mean fields if) a and ipM- The resulting predictions for the mean- field 
parameters associated with 2-atom terms are 

v M = -1/a 2 , (74a) 
qam = y/l&ir/a, (74b) 

The values of v M and qam agree with the limiting values of RDSS near the Feshbach 
resonance. Note that the value of Xaa depends on the chemical potential, which depends on 
the number density. If we take the dilute limit in which /i — > and insert the parameters 
in Eq. (!74|) into the right side of Eq. (ITT]) , we get the correct scattering length a. Note that 
the scattering length is obtained as the sum of two terms: a = (—a) + 2a. 

A comparison of Eqs. (!70|) and (!73|) shows that the truncated 1PI approximation predicts 
Aam = and Xmm = 0. These predictions are artifacts of our truncation of the 1PI effective 
action after the 2-atom terms. If the 1PI effective actions were truncated after the 3-atom 
terms, there would be an additional term XamD* A* AD in Eq. (1731) . If it were truncated 
after the 4-atom terms, there would also be a term ^Xmm(D*D) 2 . Using the exact result 
for the atom-dimer scattering length a ad in Eq. (14"5"|) . we can deduce the value of Xam- The 
T-matrix element iT for atom-dimer scattering is the sum of tree diagrams whose vertices 
are 1PI amplitudes. Thus it is the sum of the diagram in Fig. [6] whose contribution to 
the T-matrix element is given in Eq. PU|) and a diagram with a 1PI atom-dimer scattering 
amplitude. For zero collision energy, the 1PI atom-dimer scattering diagram reduces to 
—iXam and the T-matrix element reduces to T = —3tt a ad- Thus the atom-dimer scattering 
length is 

Xam 16 . , 

= ~37r y a - (75) 
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Comparing with the exact result for a ad in Eq. (|45l) . we obtain the exact result for Xam- 

\ AM = 3^(6.79- 2.15 tan [s ln(a«») + 1.16]) a, (76) 

where so = 1.00624 and the Efimov parameter. 

The location of the quantum phase transition in the mean-field approximation was de- 
termined in RPW and RDSS [l(Jlll|]. The critical number density n c at the quantum phase 
transition satisfies 1 

vm = -gAMV2n~c + (Aam - \^mm) n c . (77) 

If the dependence of the mean-field parameters on the detuning parameter, or equivalently, 
on the scattering length a, is known, this condition can be translated into a critical value 
a c of the scattering length. To compare the condition in Eq. (177)) with the results from the 
truncated 1PI effective action, we must set Xam = ^mm = 0. Using the values of Um and 
qam in Eq. (1741) . we recover the expression for the critical number density n c in Eq. (1531) . 



V. QUASIPARTICLES 

In this section, we use the truncated 1PI approximation to study the quasiparticles as- 
sociated with small fluctuations around the static homogeneous states. There are two types 
of quasiparticles: atom quasiparticles with dispersion relation u)A{k) and pair quasiparticles 
with dispersion relation ujp(k). The two types of quasiparticles can be distinguished by 
whether the asymptotic behavior of their dispersion relations as k — > oo is that for an atom 
or a pair of atoms: 

cu A (k) — >\k\ (78a) 
uj P {k) — ► \k 2 . (78b) 

If the dispersion relation uj(k) for a quasiparticle has an imaginary part, it indicates a 
dynamical instability of the system. 



A. Atom condensate 

To find the dispersion relation u(k) for quasiparticles associated with fluctuations of 
A(r, t) around the mean-field value A, we look for solutions to Eq. (130]) of the form 

A(r, t)=A + i +e -^+ ik - r + i_e + ^*- ik - r , (79) 

with A + and A_ infinitesimal. There are nontrivial solutions for A + and v4_ only if the 
dispersion relations ui{k) satisfy 

= ( G 2 c) 2 + r 2 (r-l)[/ + -/_]a 2 c-ir 2 (r-l)[/ + + /_](r 2 + a 2 A; 2 ) 

-r 4 (r - l) 2 /+/_ - ia 2 A; 2 (2r 2 + a 2 k 2 ), (80) 

where 

^ 1 — r(l + ja 2 k 2 /r 2 ± a?uj/r 2 — it) 1 / 2 ^ ^ 



1 In RPW 0, the gAM term in this equation has the wrong sign. 
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Since Eq. (1BU1) is even in u, the solutions come in pairs: if u(k) is a solution, so is —u(k). 
The invariance of T[A] under the U(l) phase symmetry A — > e JC M guarantees that there 
must be a solution that satisfies u(k = 0) = 0. This can be confirmed by observing that 
uj = k = is indeed a solution to Eq. (1801) . 

Using algebraic manipulations, Eq. (j80p can be transformed into a 6 th order polynomial 
equation in uj 2 [26| . The polynomial equation is extremely complicated, with 266 terms when 
fully expanded into a polynomial in uj, k and r. The explicit form is not very illuminating, 
so we do not write it down here. It is important that Eq. flHUj) can be transformed into a 
polynomial equation, because it guarantees that we can find all the dispersion relations that 
satisfy Eq. flHDl) . A 6th order polynomial equation in uj 2 has precisely 6 roots in the complex 
uj 2 plane. Given numerical values of r and k, the 6 roots of the polynomial equation for uj 2 
can be easily found numerically. One can then substitute the 12 roots for uj into Eq. flHUj) 
to see which ones correspond to quasiparticle dispersion relations uj(k). 

We proceed to consider various limits in which the quasiparticle dispersion relations that 
satisfy Eq. (|80|) can be obtained in closed form. 

1. Non-interacting limit 

In the limit r — > 0~, the quasiparticle dispersion relations reduce to 

u A (k) = \k 2 , (82a) 
uj p (k) = -\/ a 2 + \k 2 . (82b) 

These are the dispersion relation for an isolated atom and for an isolated dimer with binding 
energy Ed = 1/a 2 , respectively. The atom and pair dispersion relations in Eqs. ( 182]) are 
illustrated in Fig. M as dashed lines. Since they are dispersion relations for isolated atoms 
and dimers, we refer to the limit r — > 0~ as the non-interacting limit. Since uj A (k) and ujp(k) 
are real for all k, the homogeneous state is dynamically stable in the non-interacting limit. 

2. Critical density n c 

The limit r — > 1 + corresponds to the critical density n c = 1/ (327ra 3 ) defined in Eq. (1531) . 
In this limit, the quasiparticle dispersion relations reduce to 

uj A (k) = \k{k 2 + 2/a 2 ) 1/2 , (83a) 
ujp(k) = \k 2 . (83b) 

The atom dispersion relation is that for the Bogoliubov mode in a superfluid with coherence 
length a/y/2. The pair dispersion relation is that for a dimer whose binding energy Ed = 
1/a 2 is cancelled by its mean field energy. The atom and pair dispersion relations in Eqs. (|83p 
are illustrated in Fig. [9] as solid lines. Since u> A (k) and up(k) are real for all k > 0, the 
homogeneous state is dynamically stable at the critical density. 

3. Unitary limit 

The unitary limit is a — » ±oo with fixed number density n. It corresponds to r — >• ±oo 
with r/a = y/—2jl fixed. In this limit, it is convenient to introduce dimensionless scaling 
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FIG. 9: Atom and pair dispersion relations in the non-interacting limit (r — > , dashed line) and 
the critical density (r — > 1 + , solid line). 

variables for the frequency and wavenumber: a) = ua 2 /r 2 , k = ka/r. After multiplying 
Eq. flHUj) by r -4 and taking the limit r — > ±00, it reduces to 

= u 2 + (/+ - /> - + + k 2 ) - f + f- - \k 2 (2 + k 2 ), (84) 

where 

f± = 7^—^- • (85) 

Using algebraic manipulations, Eq. (184^) can be transformed into a 6 th order polynomial 
equation in Co 2 [26]. The solutions must be inserted into Eq. ( 1541 to find which ones are 
quasiparticle dispersion relations. The atom and pair dispersion relations at large k have 
expansions in powers of 1/k. 

uj A {k) = h 2 + \±i~ + ..., (86a) 
* * k 
1 - 16 

Cbpik) = -k 2 + 1 — 4- + . . . (86b) 
4 fc 4 

The sign of the imaginary part of u)A{k) is determined by the choice of branch cut for the 
square root in Eq. (1851) . The terms of order k° in Eqs. ( I86p imply that the atom and dimer 
have mean-field energies \n\ and 2|/x|, respectively. The limiting behaviors of the dispersion 
relations as k — > are 

w^(jfc) — > (0.12 ± i 0.45) A£, (87a) 
a;p(/c) — . -|(f7rn) 2 / 3 A; 2 , (87b) 
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FIG. 10: The square of the pair dispersion relation (solid line) and the real and imaginary parts of 
the square of the atom dispersion relation (dashed line and dash-dotted line) in the unitary limit 
a — ► ±00 with n fixed. 

where kp = (Q^n) 1 ^ is the Fermi wavenumber. 

The atom and pair dispersion relations in the unitary limit are illustrated in Fig. [10] by 
plotting Reju^y, Im[o)^], and Cj 2 p as functions of k. The atom dispersion relation Cja is 
complex valued for all k. This indicates the system is unstable under fluctuations of any 
wavelength. The pair dispersion relation up remains real until k decreases below the critical 
value k c « 0.907 kp, after which Cjp becomes pure imaginary. 



4- General case 

Eq. (I80il for the quasiparticle dispersion relations cannot be solved analytically for general 
values of r in the range —00 < r < and 1 < r < +00. However, analytic expression can 
be obtained for limiting values of k. At large k, the atom and pair dispersion relations have 
expansions in powers of 1/k: 

1,9 r 2 2r 2 (r — 1) . 
Mk) = \k 2 + — 2 ±i ka s - + •••> ( 88a ) 

Mk) = \k 2 + r ^- * r %~ l) +■■■■ (88b) 

The terms of order k° imply that the atom and pair quasiparticles have mean-field energies 
\r 2 Ep> and r 2 Ep), respectively. Except at the limiting values r — > 0~ and r — > 1 + , the atom 
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dispersion relation u>A{k) is complex for all k. Thus the system is unstable under nonhomo- 
geneous fluctuations with any wavelength. The pair dispersion relation up(k) remains real 
until k decreases below a critical value k c that satisfies [26i] 

= (k c af + A(k c a)\r 2 - 1) - lQ(k c a) 2 r 2 (r - 1) - 16r 4 (r - l) 2 , (89) 

and then it becomes pure imaginary. As r increases from 1 to +oo and then increases from 
— oo to 0, k c increases from to 0.907 /cf and then decreases to 0. 



B. Mixture of atom and dimer condensates 



To find the dispersion relations u(k) for quasiparticles associated with fluctuations of 
A(r, t) and D(r, t) around the mean-field values A and D, we look for solutions to Eqs. (|35|) 
of the form 

A(r,t) = A + A + e- lujt+ik - r + A_e +l ^ ik - r , (90a) 
D(r,t) = D + D + e- iuJt+ik - r + D_e +llJt - ikr , (90b) 

with A + , D + , and Z)_ infinitesimal. Surprisingly, the condition for non-trivial solutions 
can be reduced to Eq. flHUj) [26]. Thus the dispersion relations for fluctuations of A(r,t) and 
D(r, t) predicted by Eq. fl35|) are identical to the dispersion relations for fluctuations of A(r, t) 
predicted by Eq. fl30l) . This may seem quite remarkable at first, but it is a consequence of 
the equivalence of the effective actions Ti+^A; jjl] and Fi +2 [A D; /jl]. Although these effective 
actions give the same dispersion relations, the coefficients of A + , A_, D + , and D_ obtained 
by solving Eq. (131)1) provide more detailed information about the nature of the quasiparticles 
than the coefficients of A + and A_ obtained by solving Eq. (|30p . 



C. Pure dimer condensate 

If a > and < n < n c , the static homogeneous state is a dimer superfluid with mean 
fields D ^ and A — 0. The equations for small oscillations of A(r, t) around and D(r, t) 
around D decouple. The dispersion relations are 

UA (k) = [(l/(2a 2 ) + \k 2 f -8nn/a\ 1/2 , (91a) 
Up (k) = \k 2 . (91b) 

The pair dispersion relation in Eq. f!91bl) corresponds to a dimer whose binding energy 
Ed = 1/a 2 is exactly canceled by its mean-field energy. The atom dispersion relation in 
Eq. (191 al) corresponds to a mean-field energy [l/(2a 2 ) 2 — Sun/ a] 1 / 2 . The mean-field energy 
vanishes at the critical density n c = (327ra 3 ) -1 . 

The dispersion relations in Eq. (19~T|) are real-valued for all k. Thus if a > 0, the truncated 
1PI approximation predicts that the pure dimer condensate is stable for n < n c . The 
existence of Efimov states suggests that if 3-body effects were taken into account, the static 
homogeneous state containing a dimer condensate but no atom condensate would also be 
unstable for n < n c . Thus the stability of the dimer superfluid is probably an artifact of the 
truncated 1PI approximation. 
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D. Comparison with previous work 



Several groups have studied the quasiparticle dispersion relations for the many-b ody 
system of bosonic atoms near a Feshbach resonance using the mean-field approximation [lCj, 
11 . HI]. RPW and RDSS derived the quasiparticle dispersion relations in the molecular 
superfluid phase [l(J U|. Their dispersion relations reduce to our results in Eqs. (j9Tj) if we 
set Xam — ^mm = and insert the values of the mean-field parameters for the truncated 
1PI approximation given in Eqs. (IT%|) . Basu and Mueller studied the quasiparticle dispersion 



relations in the atom superfluid phase [13j. They showed that if the atom- molecule and 
molecule-molecule coupling constants Xam and Xmm are set to zero, there is a dispersion 
relation with the behavior uj 2 — *> c 2 k 2 as k — > with c 2 < 0. Thus the atom superfluid phase 
is unstable with respect to long-wavelength fluctuations. Our result for the pair dispersion 
relation up(k) in the unitary limit in Eq. (187b j) is consistent with their result. Our results 
go further by showing that the atom dispersion relation u)A{k) is complex for all k. Thus 
the atom superfluid phase is unstable with respect to fluctuations of any wavelength. 



VI. SUMMARY 

Previous work on the dimer superfluid phase of the strongly-interacting Bose gas has relied 
on the existence of a quantum field whose expectation value can indicate the presence of a 
dimer condensate. The atom superfluid phase and the dimer superfluid phase can both be 
described using the mean-field approximation. However a dimer condensate is also possible 
in models when there is no such quantum field. In this case, the mean-field approximation 
cannot be used to describe the dimer superfluid. We have developed a method for describing 
the dimer condensate in such case. It is based on the 1PI effective action T[A], which is a 
functional of classical atom field A(r, t) that encodes complete quantum information of the 
system. The solutions to the variational equation for T[A] correspond to states that contains 
an atom condensate only or a mixture of an atom condensate and a dimer condensate. The 
functional T[A] also contains information about states that contains a dimer condensate but 
no atom condensate, but it is encoded in a more subtle way through poles in functions of 
differential operators. In the case a > 0, we constructed an equivalent 1PI effective action 
T[A, D] that is also a functional of a classical dimer field D(r,t). Some of the solutions to 
the variational equation for T[A, D] correspond to the same states as the solutions to the 
variational equation for T[A}. These states contain an atom condensate only or a mixture of 
an atom condensate and a dimer condensate. However the variational equations for T[A, D] 
has additional solutions that correspond to states containing a dimer condensate but no 
atom condensate. 

To illustrate our method for describing the dimer condensate, we introduced an approx- 
imation in which 1-atom and 2-atom terms are treated exactly while irreducible iV-atom 
terms with N > 3 are ignored. This approximation can be defined by truncated 1PI ef- 
fective actions ri +2 [A] and ri +2 [A,Z)]. We determined these functionals exactly for the 
zero-range model in which the only interaction parameter is the large scattering length a. 
We demonstrated the limitations of the truncated 1PI approximation by using it to calcu- 
late the cross section for atom-dimer elastic scattering. We compared the results with those 
from exact solutions of the 3-body problem, which depend on the Efimov parameter The 
errors can be very large, so the truncated 1PI approximation is not useful as a quantitative 
approximation. However many results can be obtained analytically in this approximation, so 
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it may be useful as a benchmark against which results from other methods can be compared. 
We have found the truncated 1PI approximation to be particularly useful for illustrating 
applications of the 1PI effective action. 

We applied the truncated 1PI approximation to the homogeneous strongly-interacting 
Bose gas at zero temperature. We first determined the static homogeneous states of the 
system. The phase diagram is qualitatively consistent with that derived in Refs. [l(J 11 



There is an atom superfluid phase and a dimer superfluid phase separated by quantum 
phase transitions at a = and a c = (32rrn)~ 1 ^ 3 . The unitary limit a — > ±00 is consistent 
with Ho's universality hypothesis (25[. There is a perfectly smooth transition from an atom 
superfluid state containing a mixture of an atom condensate and a dimer condensate as 
a — ► +00 to an atom superfluid containing an atom condensate only as a —>■ —00. We also 
determined the quasiparticles dispersion relations associated with small fluctuations around 
the homogeneous states. In the atom superfluid state, the pair dispersion relation up(k) 
is pure imaginary for sufficiently small k, indicating an instability with respect to long- 
wavelength fluctuations. The atom dispersion relation u^k) is complex for all k, indicating 
an instability with respect to fluctuations of any wavelength. In the dimer superfluid state, 
the atom and pair dispersion relations are real-valued for all k, indicating that the system 
is dynamically stable. We suggested that this stability is an artifact of the truncated 1PI 
approximation. In the case a > 0, the properties of the atom superfluid phase could be 
derived either from the variational equation for or from the variational equations for 

ri +2 L4, D]. They give exactly the same thermodynamic properties and identical quasiparti- 
cle dispersion relations. This illustrates the equivalence of these two truncated 1PI effective 
actions. Some advantages of ri +2 [A, D] are that it allows the non-condensate contributions 
to thermodynamic properties to be separated into a dimer-condensate contribution and a 
remainder and it gives more detailed information about the nature of the quasiparticles. The 
properties of the dimer superfluid phase can be derived easily from the variational equations 
for ri + 2[v4, D]. That information must also be encoded in ri+2[A], since these two functional 
are equivalent, but is is in a much less accessible form. 

Much of the previous work on the Bose gas near a Feshbach resonance has been carried 
out with the mean- field approximation using the Hamiltonian density in Eq. (|70|) . which 
depends on atomic and molecular mean fields ipA and ipM- We used the truncated 1PI 
effective action r 1+2 [v4, D] to deduce the limiting values of some of the mean-field parameters 
when the system is sufficiently close enough to the Feshbach resonance. Our results for 
and gAM agree with those of Ref. [ll|. Our result for Xaa increases linearly with a, while 



Xaa was assumed to be constant in Ref. 111 . We also used results from exact solutions to 



the 3-body problem to deduce the limiting value of Xam- Our result for A am in Eq. (ITS]) 
has log-periodic dependence on the Efimov parameter. The simple result Xam —* 327ra in 
Ref. [II] was derived using an approximation equivalent to the truncated 1PI approximation. 

The method we have used to describe states containing a dimer condensate can be gen- 
eralized to states containing condensates of Efimov trimers. Information about a specific 
Efimov trimer is encoded in the 1PI effective action V[A] through the poles of functions of 
differential operators that act on products of three classical atom fields. By introducing a 
classical trimer field T(r,t), it should be possible to construct an equivalent 1PI effective 
action T[A, T] or T[A, D, T] in which these poles have been eliminated. One could then use 
the variational equation for these effective actions to study the trimer superfluid phase of 
the strongly-interacting Bose gas. 
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